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We study the effect of an external magnetic field in the Kondo regime of a double-quantum-dot 
system in which a strongly correlated dot (the "hanging dot") is coupled to a second, noninteracting 
dot that also bridges the gap between two external leads. In zero field, the spectral function of the 
hanging dot has previously been shown to exhibit a split-peak structure near the Fermi level due 
to "Kondo resonance filtering" by the bridging dot. We show using the numerical renormalization 
group that application of a magnetic field leads to a subtle interplay between electronic interference, 
Kondo physics, and Zeeman splitting with nontrivial consequences for the spectral and transport 
properties. The value of the hanging-dot spectral function at the Fermi level exhibits a nonuniversal 
field dependence that can be explained using a generalized Friedel sum rule for a Kondo system 
with energy-dependent hybridization. The magnetic field also accentuates the exchange-mediated 
interdot coupling, which dominates the ground state at intermediate fields leading to the formation 
of antiparallel magnetic moments on the dots. By tuning gate voltages and the magnetic field, 
one can achieve complete spin polarization of the linear conductance between the leads, raising 
the prospect of applications of the device as a highly tunable spin filter. The system's low-energy 
properties are qualitatively unchanged by the presence of weak on-site Coulomb repulsion within 
the bridging dot. 

PACS numbers: 73.63.Kv, 72.10.Fk, 72.15.Qm 



I. INTRODUCTION 



Electron correlations in quantum-dot structures re- 
sult in many fascinating effects that can be probed 
in detail with remarkable experimental control of sys- 
tem parameters Perhaps one of the most interesting 
regimes occurs when electrons confined in the dot acquire 
antifcrromagnetic correlations with electrons in the leads, 
giving rise to the well-known Kondo effect^ The sim- 
plest realization of this phenomenon in a single quantum 
dot is characterized by just one low-energy scale, set by 
the Kondo temperature, that controls (among other fea- 
tures) the width of a many-body resonance at the Fermi 
energy^ Recent experimental 5 - - — and theoretical^— 
studies of the Kondo effect in multiple quantum dots 
have revealed a complex competition between geometry 
and correlations, making evident that these structures 
provide a flexible setting in which to explore much novel 
physics. 

In this context, we have shown that certain double- 
quantum-dot devices exhibit unexpected manifestations 
of Kondo physics, and have predicted conductance signa- 
tures of these effects in realistic experimental setups i 20 ' 21 
A particularly striking situation arises when a small, 
strongly interacting "dot 1" is connected to external leads 
via a large "dot 2" that is tuned to have a single-particle 
level in resonance with the common Fermi energy of the 



leads. In this "side-coupled" or "hanging-dot" configu- 
ration, the Kondo resonance, which normally has a sin- 
gle peak at the Fermi energy, splits into two peaks — a 
behavior that can be understood as a consequence of in- 
terference between the many-body Kondo state in dot 1 
and a singlc-particlc-like resonance that controls (or "fil- 
ters") its connection to the leads.— ~— The magnitude 
of the Kondo peak splitting is determined by the bal- 
ance of several important energy scales in the problem: 
the width and position of the active single-particle level 
in the dot 2 that acts as a bridge between the leads; the 
height of the effective single-particle resonance set by the 
interdot coupling; and the many-body Kondo tempera- 
ture, which is determined by the preceding energy scales 
in combination with the dot-1 level-position and interac- 
tion strength. This filtering of the leads preserves a fully 
screened Kondo ground state with a Kondo temperature 
that rises with increasing interdot coupling. 

In this work, we investigate the effects of an 
external in-plane magnetic field on such a double- 
quantum dot-system in the hanging-dot arrangement. 
The field — which introduces another energy scale, the 
Zeeman energy — is known to be detrimental to the 
Kondo state in single-dot systems Jj22r— Using numerical 
renormalization-group methodsj^ 4 ^ we study the inter- 
play between the different energy scales and discuss the 
behavior of the Kondo resonance in the presence of com- 
peting interactions. This interplay reveals itself in the 



fundamental Fermi-liquid properties of the system, such 
as the variation with magnetic field B at zero tempera- 
ture of the Fermi-energy (<j = 0) value of the hanging-dot 
spectral function A±(u},T). Instead of the usual mono- 
tonic deca y 25 ' 26 of Ai(Q, 0) with increasing B we find a 
markedly nonuniversal behavior, where Ai(0,0) passes 
through a maximum at a nonzero value of the field. This 
effective field- enhancement of the Kondo spectral func- 
tion is a consequence of the hanging-dot geometry. The 
same behavior can also be understood using a general- 
ized Friedel sum rule for a Kondo system with energy- 
dependent tunneling to/from the leads. The sum rule 
predicts parameter- and field-dependent phase shifts that 
impart the unusual nonmonotonicity to the variation of 
.Ai(0,0) with B. 

In addition, we show that the competition between 
Zeeman splitting of the dot levels and Kondo screen- 
ing results in a dominant exchange-mediated antiferro- 
magnetic coupling of the dots over a range of moderate 
magnetic fields, before both dots become fully polarized 
at higher fields. Finally, we report interesting transport 
properties that arise from the aforementioned phenom- 
ena. A key result is the generation of spin-polarized cur- 
rents through the device, which can be tuned by adjust- 
ing gate voltages to achieve total polarization. 

The remainder of the paper is organized as follows: In 
Scc.|lT]wc describe the effective Anderson impurity model 
for the double-quantum-dot system. Section IIIII presents 
the low-energy spectral properties, while Sec. IIVI inter- 
prets the nonuniversal behavior of Ai(u> = 0, T = 0) vs 
B in terms of a generalized Friedel sum rule. The trans- 
port properties, including spin polarization, are explored 
in Sec. |VJ Concluding remarks appear in Sec. IVII 
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II. DOUBLE-QUANTUM-DOT SYSTEM 



The system under study, which is depicted schemati- 
cally in Fig. [TJ contains two quantum dots. Dot 1 has 
a large Coulomb repulsion U\ when its single active en- 
ergy level is doubly occupied. Electrons in this level can 
tunnel only to/from dot 2 with a matrix element A. Dot 
2 has negligible electron-electron interactions (U2 — 0) 
and one active level that can be tuned by gate voltages 
to be at or near resonance with the common Fermi en- 
ergy €f = of left (L) and right (R) leads. In addition 
to their coupling to dot 1, electrons in dot 2 can tunnel 
to/from each lead with a hybridization matrix element 
V 2 . 

The system can be described by a variant of the two- 
impurity Anderson Hamiltonian: 
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FIG. 1: (Color online) Schematic representation of the side- 
coupled double-quantum-dot system. The dot QD1 has a 
strong Coulomb interaction Ui and is coupled only to the 
second dot, labeled QD2. The latter dot has negligible local 
interactions (i.e., U2 — 0) and the energy of its active level is 
tuned to allow tunneling at or near resonance with the Fermi 
level of the left (L) and right (R) leads. 
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Here, di„ annihilates an electron in dot i with spin z 
component ^cr (a = ±1 or cquivalently j", l) and en- 
ergy e l<J — Si + \agi[i B B ; n la = d\ a d lG is the cor- 
responding number operator; and Qkc annihilates an 
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with spin z component ^er and en- 
ergy EfkCT = ££k + ^crgcHsB. The magnetic field Bz 
with B > is assumed to lie in the plane of the two- 
dimensional electron gas in which the dots and leads are 
defined, so that it produces no kinematic effects and en- 
ters only through Zeeman level splittings. This Hamilto- 
nian differs from a generic two-impurity Anderson model 
through the absence of a dot-1 hybridization Vi, a con- 
sequence of the hanging-dot geometry. Throughout the 
greater part of the paper, we also take U2 — 0, a case that 
is particularly convenient for algebraic analysis. The ef- 
fect of nonvanishing dot-2 interactions is addressed at the 
end of Sec. El 

Without loss of generality, we take all tunneling matrix 
elements to be real. We consider local (k-independent) 
and symmetric coupling to both leads, and assume that 
the dots have equal effective g factors g\ = g 2 = g- sim- 
plifications that do not qualitatively affect the physics. 
The leads are assumed to have featureless band struc- 
tures near the Fermi energy, modeled by the flat-top den- 
sities of states = pr(w) = p{oS) = (2D)^ 1 Q(D — 
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\oj\) where D is the half-bandwidth and Q(x) is the Heav- 
iside function. In this situation, the Zeeman splitting 
of the conduction band produces only very small effects 
near the band edges, so for convenience we set the bulk 
g factor to g c = throughout what follows. 

The primary quantities of interest in this work 
arc the retarded dot Green's functions Qi a {u>,T) = 

((d Ur ;4*))» fOT * = i.2, where ((A;B)) U = 
-i f™({A(t),B(0)})e ZUJt dt and (• • • ) denotes an appro- 
priate thermal averaged In particular, we are interested 
in the spectral functions Ai a (u),T) = — 7r~ 1 Imf/.; cr (w, T) 
and the system's linear (zero-bias) conductance, given 
within the Landauer approach as G = G a with 



GAT) 



-lm%(u,T)](-df/dw)du, (5) 



where /(w,T) is the Fermi distribution function at tem- 
perature T. In the side-connected geometry, the trans- 
mission is^i 



where A 2 = nV 2 2 /D. Thus, 



(6) 



G a (T) = ^-ttA 2 f A 2a {u, T) (-df/dw) du, (7) 



which reduces at zero temperature to 

e 2 

G a (T = 0) = -nA 2 A 2a (0,0). 



(8) 



In order calculate the dot spectral functions A ia (uj 7 T) 
taking full account of the strong electronic correlations 
arising from the U\ term in Eq. ([2]), we employ the numer- 
ical renormalization-group (NRG) method, performing a 
logarithmic discretization of the states in the leads and 
iteratively solving the discretized Hamiltonian. In eval- 
uating the spectral functions, we perform a Gaussian- 
logarithmic broadening of discrete poles obtained by the 
procedure described in Ref. H3- At temperatures T > 
we use the density-matrix variant of the NRG^i which 
has better spectral resolution at high frequencies and 
nonzero fields.—^ 

All numerical results were obtained for a symmetric 
dot 1 described by Ui = -2ei, for dot 2 width A 2 = 0.02, 
and for NRG discretization parameter A = 2.5. Except 
where it is stated otherwise, we consider a strongly cor- 
related dot 1 with Ui = 0.5 and situations in which a 
nonintcracting dot 2 is tuned to be in resonance with 
the leads, i.e., U 2 = £ 2 = 0. We adopt units in which 
D = h = ks = g^B = 1- 



III. SPECTRAL PROPERTIES 

As shown previously^ the side-dot arrangement can 
lead to a strong Kondo state in the hanging dot with an 
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FIG. 2: (Color online) (a) Spin-averaged dot-1 spectral func- 
tion Ai vs frequency ui at zero temperature for JJ\ = — 2ei = 
0.5, £2 = 0, A = 0.0627, and (from bottom to top curve, 
offset for clarity) B = 0, 0.01, 0.02, and 0.03. The spectral 
function is multiplied by tiTi where I\ = A 2 /A 2 . Insets: Ex- 
panded views of A\ vs ui around the Fermi level u = for 
the same system, with B ranging from to 0.035 (bottom 
to top, curves offset for clarity) in steps of 0.005 in the left 
inset, and from 0.04 to 0.07 (bottom to top, curves offset for 
clarity) in steps of 0.01 in the right inset, (b) Spin-up (black 
squares) and spin-down (red circles) dot-1 spectral functions 
Aicr(io,T = 0) for B — 0.02 with all other parameters as in 
(a), (c) Same as (b), except for B — 0.04. 



unusual spectral function Ai a (uj,T) that in most cases 
exhibits a split-peak structure at low temperatures. This 
split Kondo resonance is, however, accompanied by the 
complete screening of the magnetic moment on dot 1, 
resulting in a Kondo singlet ground-state. The splitting 
can be understood^ , from a mapping (detailed in Scc. lIV[) 
of the Hamiltonian ([1]) onto a single-impurity Anderson 
model with an energy-dependent hybridization function 
having a peak of width A 2 and height proportional to A 2 . 
This energy dependence gives rise to an intrinsic splitting 
of the spectral function that grows with the interdot cou- 
pling. In contrast with the field-induced splitting of the 
Kondo resonance in the conventional Anderson model, 
the zero-field splitting does not signal a weakening of the 
screening by the leads. Rather, it can be interpreted as a 
consequence of interference between (i) the single-particle 
resonance of width A 2 represented by the effective hy- 
bridization function and (ii) the quasiparticle resonance 
at the Fermi level due to the Kondo effect, whose width is 
given by the Kondo temperature Tjf. The robustness of 
the Kondo effect is shown by a rise of Tk with increasing 
A. The zero-field splitting is spin-symmetric since it does 
not arise from breaking of time-reversal symmetryj^ 2 - 

The Kondo-peak splitting described in the previous 
paragraph has superficial similarities to splittings found 
in a single quantum dot. In the latter context, one can 
think of a field as displacing the up- and down-spin spec- 
tral functions by opposite amounts Aw ff = —^ag/isB. 
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Once B > Tx/g[J-B, the up- and down-spin spectral res- 
onances have little overlap, inhibiting the coherent spin- 
flip processes that are responsible for the Kondo effect. 
Thus the facts that the spin-averaged spectral function 
splits around the Fermi level and its value at the Fermi 
energy undergoes a monotonic decrease with increasing 
field signal a progressive weakening of the Kondo effect. 

Having seen that the zero-field splitting in the spectral 
function of a hanging dot and the splitting of the single- 
dot Kondo resonance in the presence of a magnetic field 
(or ferromagnetic lead s 24 ' 29 ) represent very different phe- 
nomena, we are led to consider the combination of these 
effects when a Zeeman field is applied to double quan- 
tum dots. Figure dja) shows the spin-averaged spectral 
function 

A 1 (u,T) = l[A n {u,T) + A u (u,T)] (9) 

for a hanging-dot setup at zero temperature with U\ = 
0.5 and A = 0.0627. The different curves, vertically offset 
for clarity, correspond to four different values of B. For 
zero field (the bottom curve), Aif(w,Q) = Ai±(ui,0) = 
Ai(u>,0), so each spin- resolved spectral function shows a 
symmetric Kondo-peak splitting due to the interdot cou- 
pling A. With increasing B, the split peaks merge into 
a single peak at uj = 0, clearly seen for B = 0.03 (top 
curve) . The left inset to Fig. [2ja) shows in greater detail 
the convergence of the peaks near the Fermi energy, with 
the maximum in A\ vs uj at uj = being best defined 
at B = 0.035, a field where, incidentally, the absolute 
value of Ai(0,0) exceeds that at B = by nearly a fac- 
tor of two. For slightly larger fields, the central peak 
again splits into two before all low-energy features be- 
come flattened out at fields B > 0.07 [right inset to Fig. 
Ha)]. 

The field-induced merging of the peaks in Ai(ui,0) 
arises from opposite displacements of A\^{ui, 0) and 
A\i(uj,Q) along the uj axis. In a nonzero magnetic field, 
A la (uj,T) ^ A la (-uj,T) but Ai t (w,T) = A u (-u,T). 
This is illustrated for B = 0.02 in Fig. HJb), which also 
shows that the heights of the two peaks in each spin- 
resolved spectral function Ai a (uj,0) are no longer equal. 
Upon further increase in the field to B = 0.04 [Fig.^c)], 
the double-peak structure is replaced by a single peak 
near uj = in each spin-resolved spectral function. For 
larger values of B, these peaks move away from the Fermi 
energy and the usual Zeeman-splitting of the Kondo peak 
with decreasing amplitude becomes evident in the spin- 
averaged spectral function [right inset in Fig. [Ha)]. 

We now focus on the field dependence of A\(uj = 0, T = 
0), a quantity that acts as a sensitive measure of the in- 
terplay of the different energy scales in the problem: the 
single-particle resonance width A2, the Kondo tempera- 
ture Tk, and the Zeeman energy g^sB. Figure[3{a) plots 
7rA(0) Ai(0,0) vs B/Tk (taking g^s = 1) for six values 
of A. The energy scale A(0) will be defined in the next 
section; for now, it suffices to note that it is proportional 
to [I + (B/2A2) 2 ] -1 , i.e., it is a decreasing function of the 
field. The figure reveals two distinct regimes of behavior: 




B/T K 

FIG. 3: (Color online) (a) Spin-averaged dot-1 spectral func- 
tion at the Fermi level Ai(uj = 0, T — 0) vs scaled magnetic 
field B /Tk for Ui = — 2ei = 0.5, £2 = 0, and six values 
of A. Ai(0, 0) is multiplied by the field-dependent quantity 
7rA(0) defined in Eq. (|26[) . The larger A values produce a 
nonmonotonic field variation of Ai(0, 0), with a peak around 
B ~ 2Tk- Inset: Corresponding plot for the noninteracting 
case Ui — £1 = 0, with the field scaled by the interdot cou- 
pling A. (b) Phase factor ipi^ — — tp 1 _i corresponding to the 
data in (a), determined from the generalized Friedel sum rule 
Eq. Q27p using the magnetization data plotted in Fig. [J] 



(1) For A < 0.05, 7rA(0) Ax(0, 0) decreases monotonically 
from its zero-field value 1 over a characteristic field scale 
that increases with A (and is not simply Tk, as it is in 
the single-dot case). (2) For A > 0.05, ttA(0) Ai(0, 0) 
has a striking nonmonotonic variation with increasing 
B, reaching a second maximum 7rA(0) Ai(0, 0) = 1 at 
B = B* ~ 2Tk, beyond which field it decreases mono- 
tonically. In view of the field dependence of A(0), the 
value of Ai(0, 0) at B = B* is [I + (5*/2A 2 ) 2 ] times its 
zero-field counterpart. The two regimes seen in Fig.[3ja) 
are in sharp contrast with the monotonically decreasing 
and universal dependence of the Fermi-energy spectral 
function on B/Tk in the conventional single- impurity 
Kondo^ and Anderson 2 ^ models. The next section will 
analyze these behaviors in terms of a generalized Friedel 
sum rule. 



IV. GENERALIZED FRIEDEL SUM RULE 

In Scc. lIV Al we review the Fcrmi-liquid relation known 
as the Friedel sum rule 4 - that sets the Fermi-energy value 
of the zero-temperature spectral function in the one- 
impurity Anderson model, and discuss its generalization 
to problems with nontrivial structure in the density of 
states. Section HV B I explains how this generalized Friedel 
sum rule accounts for the variation of A\(uj = 0, T = 0) 
in our doublc-quantum-dot system. 
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A. Single Anderson impurity 

In this subsection, we consider a single-impurity An- 
derson model 



H = ^ SdaTlda + U UdfUdi + ^ ^kerC^C^ 



k,cr 



E(^4c kCT + H.c), 



(10) 



k,<r 



specified in Eq. (JTTJ) , but YP dtT (ui,T) is the retarded im- 
purity self-energy for the noninteracting system [Eq. (fTU)) 
with [/ = 0], which satisfies ImS^(u,T) = -A CT (w). 

In order to allow for the most general situation where 
A-|-(w) 7^ A^(w), rather than considering the spin- 
averaged spectral function, it is convenient to focus on 
the dimensionless, hybridization-weighted spin-average 

%T) = ^A ff ( U )4( W ,r) (16) 



where e dr , = e d + \ag[i B B and e kcr = e k + \ag c ^ B B. 
The conduction-band dispersion e k and the hybridiza- 
tion Vic enter the impurity properties only in a sin- 
gle combination: the zero-field hybridization function 
A (w) = 7r^ k |Vk| 2 <$(w — £k). We denote the fully inter- 
acting retarded impurity Green's for this problem by 



g d(T (uj,T) = ((d a ;dt)) u = 



1 



'da 



E da (uj,T) ' 

(11) 

where T. da (uj,T) is the retarded impurity self-energy. 

In the conventional Anderson model, where the hy- 
bridization function is assumed to take a flat-top form 
A (w) = TQ(D — \uj\), the Friedel sum rule relates 
the Fermi-energy value of the zero-temperature, zero- 
field impurity spectral function A d (cu,0) = Ada(u,Q) = 
— Tr^^inQdai^i 0) to the average impurity occupancy 
(rid) = (nd-\) + (ridi) as 



7^(0,0) 



sm 



2 {Ud) 



(12) 



In the special case Ed = ~\U of a symmetric impu- 
rity, where particle-hole symmetry ensures that (n d ^) + 
( n di) = 1, Eq. (fT2j) has been extended^! to show that 
in a Zeeman field _B, the spin-averaged impurity spectral 
function A d (w,T) = MA d f(uJ, T) + A di (uj,T)} satisfies 



7rr^ d (0, 0) = cos 2 irM d . 



(13) 



where M d (B) — \((n d f) — (ndi)) 1S ^ e impurity magne- 
tization in units of gfi B . 

Our goal is to generalize Eqs. (|T2"j) and ([TB")) to al- 
low for arbitrary impurity parameters and any form of 
Ao(w). One can sho w 4 i 22 i 31 that provided the system 
is in a Fermi-liquid regime [where the imaginary part of 
T, da (uj,T = 0) varies as uS 2 for ui —5- 0], the spin-resolved 
spectral functions at zero temperature satisfy 



ttA ct (0) A d(7 (0,0) 



! (vr(n d(T ) + v3 CT ), (14) 



where A a (ui) = Aq(uj — ^ag c /j, B B) and 
'° 5£°>,T = 0) 



fa = Im 



duj 



g d(T (( J j,T = 0)du> (15) 



is a spin-dependent phase shift. In Eq. (fT51 . Q dcr (u>,T) is 
the fully interacting retarded impurity Green's function 



with the useful property that the (spin-summed) linear 
conductance through a single quantum dot described by 
the Anderson impurity model is 



G(T) 



2e*_ 
h 



F(uj,T)(-df/duj)duj. (17) 



with a zero-temperature limit 



G(T = 0) 



2e 2 



F(0,0). 



(18) 



The hybridization-weighted function reduces to a mul- 
tiple F(ui,T) = 7rA(w) Ad(u), T) of the spin-averaged 
spectral function if A^(w) = Aj_(u;) = A(w), a condi- 
tion that is satisfied (i) for all values of uj in zero mag- 
netic field, and (ii) at ui = for any field B such that 
A (\g c ^ B B) = A (~y ct i B B). 

Inserting Eq. ([H]) into Eq. (fTI)]) , rewriting (nda) = 
\(nd) + &M d , and defining tp± = ip^± ip^, one obtains 

F(0,0) = ![l-cos(7r(n d )+¥>+) cos(27rA/ d + ^_)] . (19) 

This generalized Friedel sum rule relates the value of 
the hybridization-weighted spin-averaged spectral func- 
tion at oj = and T = to the impurity occupancy, 
the impurity magnetization, and spin-dependent phase 
factors that account for the energy dependence of the 
hybridization function. The right-hand side of Eq. (|T9]) 
has a maximum possible value of 1 , implying through Eq. 
([T5]l that G(T — 0) < 2e 2 //i, as one would expect for a 
problem with a single transmission mode in the left and 
right leads. 

In general, each of the phase factors tp^ and ip^ has 
a complicated dependence on the hybridization function 
A (a;), the impurity parameters U and e d , and the mag- 
netic field B. This makes it highly improbable that for a 
generic choice of model parameters there exists any value 
of B for which the system satisfies the requirements 

cos(7r(n d ) +ip + ) = - cos(27rAf d + p_) = ±1 (20) 

for achieving F(0, 0) = 1 and, hence, a unitary conduc- 
tance G(T = 0) = 2e 2 /h. 

However, under conditions of strict particle-hole sym- 
metry, where £ d = ~\U and each conduction-band state 
k having e k > can be put into one-to-one correspon- 
dence with a state k(k) having e k ( k ) = — £ k , the Hamilto- 

nian |T]) is invariant under a transformation d a — > ±d_ a , 
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c kcr — > T C £._ CT - This invariance leads to the symme- 
tries A t (w)'= A^-w), E^w.T) = -^(-w.T)]* 
and Cy^(w,T) = — [<7dj,(— w, which in turn im- 

ply that A dt (oj,T) = A dir (-u),T) and ^ = -tp± (or 
<£>+ = 0). Since the particle- hole symmetry also ensures 
that Aq(w) = Ao(— ui) and (n^) = 1. it follows that 
A-f-(O) = A 4.(0) and that the generalized Fricdel sum rule 
reduces to 



7rA(0)A d (0,0) 



cos 



(■7rM d + tpf). 



(21) 



In situations described by Eq. ([21) , the conductance will 
reach its unitary limit whenever (nMd + Vt)/ 71 " equals an 
integer. It is much more likely that this single condition 
will be met under variation of B than that the system 
will simultaneously satisfy both parts of the condition 
Eq. (|20[) applicable away from particle-hole symmetry. 

The conventional fiat-top hybridization function 
Ao(w) = T<d(D — \u>\) is not only particle- hole symmet- 
ric, but — in the wide-band limit where D greatly exceeds 
all other energy scales in the problem — yields vanishingly 
small values of tp a , thereby simplifying Eq. (|19[) to 



nrA d (0, 0) = i [1 - cos(n(n d )) coa(2irM d )] . (22) 

If the impurity is also at particle- hole symmetry (i.e., 
Ed = —hll), then substitution of <p^ — into Eq. 
(f2"Tj) . or of (rid) = 1 into Eq. (|2"2"j) . recovers the previ- 
ously derived Eq. (fTSl . Even for Ed ^ \U ■, one expects 
\Md(B) | to be an increasing function with a limiting value 
\Md(B — > oo)| = |, and therefore Ad(0, 0) to decrease 
monotonically with increasing B. 



B. Double quantum dots 

We now return to the double-quantum-dot setup de- 
fined in Eq. (UJ). It has been shown2S that for the special 
case U 2 = of a nonintcracting dot 2, the properties of 
dot 1 arc identical to those of the impurity in a single- 
impurity Anderson model of the form of Eq. (|10[) with 
U = Ui, Ed = £1, and a zero- field hybridization function 



where 



Ao(w) = n\ 2 p 2 (u;), 



1 



p 2 (ui) = - 



7T (UJ 



: 2 ) 2 + A2 



(23) 



(24) 



describes a unit-normalized Lorentzian resonance of 
width A2 [defined after Eq. ©] centered on energy 
ui = e 2 - 

In a Zeeman field B, the spin-dependent hybridization 
function of the effective one-impurity problem is 

Acr(w) = A (w — \og\i B B)- (25) 

For the resonant case £2 = considered in Figs. [5] and [31 

A (0) 




FIG. 4: (Color online) Magnetization of dot 1 (empty sym- 
bols) and dot 2 (filled symbols) vs scaled magnetic field B/Tk 
at zero temperature for the same parameters as in the main 
panels of Fig. [3] The dot-1 magnetization M\ decreases mono- 
tonically from zero over a characteristic field scale that grows 
with A and approaches 2Tk for sufficiently large interdot cou- 
plings. The dot-2 magnetization A/2 is of opposite sign to 
Mi for B < 2Tk , pointing to the dominance of the antiferro- 
magnetic interdot exchange interaction over this field range. 
Both dots become fully polarized antiparallel to the field for 
B 3> ITk- Inset: M2 vs B/X for the noninteracting system 
with the same parameters as in the inset of Fig. [3ja). In 
contrast to the interacting case, Mi decreases monotonically 
from zero with increasing field. 



in units where g^B = F Taking into account also 
the particle-hole symmetry present for E\ = ~\U\ and 
£2 = 0, the generalized Friedel sum rule [Eq. (f2Tj) ] gives 
(after translation back into the variables of the double- 
dot problem) 

/!(5) = ^M0)>li(0,0) 



2 

cos 2 (7rMi + tpi-\), 



(27) 



A t (0) = A ; (0) = A(0) 



1 + (B/2A 2 y 



(26) 



where Mj = §((?ti-|-) — (^4.)) is the magnetic moment on 
dot i, and 

J-00 ouj 

with lmY,° 1(T (u>,T) = -A a (uj). 

FigurcUJshows the variation of Mi with B for the same 
model parameters as were used in Fig. [5J As expected, 
Mi decreases monotonically from zero over a field scale 
that grows with A. For small A, this scale is identical 
to that characterizing the initial decrease of f± from 1 
[see Fig. [3](a)], while for larger A, \Mi \ grows on the scale 
B* of the second peak in fi{B). In all cases, dot 1 is 
essentially fully polarized for B > 2Tk- That the mono- 
tonic evolution of Mi does not accompany a monotonic 
decrease in fi(B) is an indication of the importance of 
the phase factor cpi^ on the right-hand side of Eq. (|27p. 
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It is difficult to evaluate tp\ a directly from Eq. (|T5l) us- 
ing the NRG because this task requires accurate determi- 
nation of both the real and imaginary parts of Q a (u>,0) 
for all u) < 0, whereas the NRG is well-suited only to 
computation of ImQ a (ui,0) for <C D. At particle- 
hole symmetry, however, one can use Eq. (|27[) to work 
backward from the NRG values of /i and Mi to find 
tfi^ = —<fi\,. Fig. (3Jb) plots the phase obtained in this 
manner from the data in Figs.[3ja) and|4j For all values 
of A, ip\^ is zero at B — (as expected) and approaches 
7r at large field values. For larger values of A, ipi^ shows 
a pronounced kink at B = B* . This kink is related, via 
Eq. (|27l), to the peak in f x (B) at B* , since M X {B) is a 
smooth function of B (as shown in Fig. |4]). 

Figure |4] also plots the field dependence of the dot-2 
magnetization. The fact that M 2 is of opposite sign to 
Mi for B < 2Tk indicates that the interactions in dot 1 
combine with the interdot hopping to yield a dominant 
antifcrromagnetic interdot exchange interaction. Over 
this range of B, it appears that the system minimizes 
its energy by first aligning the partially Kondo-screened 
magnetic moment of the strongly interacting dot 1 along 
the direction favored by the field, and then orienting the 
less-developed moment on dot 2 to minimize the inter- 
dot exchange energy even at a cost in Zeeman energy. 
The data show that this tendency becomes weaker for 
stronger interdot couplings, presumably because the in- 
terdot exchange ~ A 2 grows more slowly than the energy 
scale Tk for breaking the Kondo singlet. For all values of 
A, once B > 2Tk, the Zeeman field has largely destroyed 
the Kondo effect, and both dots are fully polarized for 
B > 2T K . 

One can gain further insight into the results presented 
in Figs. [3] and |4] by considering the limit where both dots 
are noninteracting. Equations ([TO]) and (|2"T)) hold equally 
well for interacting and noninteracting problems. How- 
ever, the case U\ = U2 = offers the advantage that 
Ai(0, 0) can also be calculated directly from the imagi- 
nary part of 



1 



u + i0+ -e lCT -E°(o;,T)' 



(29) 



where at zero temperature the noninteracting self-energy 
is 



E>, OH [(w - e 2ff )/A 2 - z]A ff ( w ), 



giving 



1 



Mo) 



7T [e 2(T A (T (0)/A2 -ei CT ] 2 + [A CT (0)] 5 



(30) 



(31) 



The hybridization- weighted spin-average of ^^(0, 0) sat- 
isfies 



h[B) -\^l + (e 2 l-e la ) 



(32) 



where e lcr = ei cr /A cr (0) and e 2cr = £2<r/A 2 . It should be 
noted that e\ a depends on B both through the Zeeman 



shift of £\ and the value of A CT (0) = Ao(— \oB). From 
Eq. (|32[) it is apparent that fi{B) attains its maximum 
value of 1 only if e 2cr = e\ a for both spin orientations, a 
condition that can be satisfied only for e± = £2 = and 
either B = or (if A > A 2 ) B = B* = 2^/A 2 - A 2 . For 
£i/0 and/or e 2 ^ 0, fi(B) may have zero, one or two 
maxima at nonzero fields, but f\ < 1 for all B. These ob- 
servations are consistent with the conclusion drawn from 
the generalized Friedel sum rule that f± = 1 is likely to 
be achieved only under conditions of strict particle-hole 
symmetry. 

The inset of Fig. EJa) illustrates the field variation of 
fi for the particle-hole-symmetric case U\ = E\ = e 2 = 0, 
with all other parameters as in the main panel. For each 
of the A values illustrated (all of which lie in the range 
A > A 2 ), fi reaches 1 at a magnetic field consistent with 
the value B* derived in the previous paragraph. Note 
that B* approaches 2A from below in the limit of strong 
interdot coupling. The inset of Fig. |4] plots M 2 vs B for 
the same noninteracting cases. For each A value, |M 2 | 
shows a purely monotonic field variation, with a rather 
sudden increase around B ~ 2A, a behavior that is mim- 
icked in the interacting system for B ~ 2Tk , especially 
at large interdot coupling A. The variation of the inter- 
acting /1 and M 2 for B > 2Tk seen in the main panels 
of Figs.[3Ja) and 2J particularly for the larger values of A, 
may perhaps be interpreted as a many-body analog of the 
noninteracting behavior in the insets, with Tk serving as 
a renormalizcd value of the single-particle scale A. 



V. ELECTRICAL CONDUCTANCE 

While the spectral functions discussed in the preced- 
ing sections are difficult to access directly in experiments, 
they may be probed indirectly through transport mea- 
surements. In this section, we show that even though the 
interacting dot 1 is not directly connected to the external 
leads, the zero-bias electrical conductance through the 
double-dot device contains clear signatures of the nonuni- 
versal variation of 7rA(0) A\(0, 0) with applied field. In 
particular, we demonstrate the feasibility of generating 
currents through the system that are strongly or even 
completely spin-polarized. 

Although the linear conductance is given most com- 
pactly by Eq. (jTJ) , it is also useful to express G in terms 
of the Green's function for the interacting dot 1 by com- 
bining Eq. ([5]) with a generalization of Eq. (6) in Ref. [2l] 
to include the Zeeman field: 

-Im7^(w,T) 

= [1 - 27rA 2 p 2(T (w)] nA a (u) A la (cj, T) + 7rA 2 p 2(T (w) 
+ 2ir(u - E20) P 2*{u) M") Re£i CT (w,T), (33) 

where p 2cr (w) = J o 2 (w— ^ag^sB), with p 2 ( w ) an d A„(w) 
as defined in Eqs. (f!H]) and (|25p . respectively. The term 
7rA 2j o 2o -(oj) describes the bare transmission through dot 2 



FIG. 5: (Color online) Linear conductance G vs scaled mag- 
netic field B /Tk at zero temperature for the same parameters 
as in the main panels of Fig. G rises from zero over the 
same characteristic field scale as governs the rise of |Mi| in 
the main panel of Fig. [3] 



in the absence of dot 1, while the remaining terms repre- 
sent additional contributions arising from conductance 
paths that include dot 1. In the special case A = 
where the latter contributions necessarily vanish, the 
zero-temperature conductance reduces to 

G onc . dot (T = 0) = C - £ — V , (34) 
where e 2(T is as defined after Eq. (|32|). 



A. Zero temperature 

Figurc[5]plots the zero-temperature linear conductance 
G as a function of scaled field B/Tk for the same param- 
eters as used in Fig. [3] For the case £2 = considered 
here, the conductance of dot 2 alone, G nc-dot(T = 0) = 
(2e 2 /h)[l + (B/2A2) 2 ] -1 , decreases monotonically from 
its unitary limit as the Zecman field detunes the dot level 
from the Fermi energy of the leads. For any A 7^ 0, by 
contrast, it has been shown previously^ that at zero field 
the presence of Kondo correlations in dot 1 leads to per- 
fect destructive interference with the resonant conduc- 
tance channel through dot 2 [i.e., complete cancellation 
at L) = of 7rA2/?2tT by the other terms on the right-hand 
side of Eq. (|33[) ]. resulting in zero conductance through 
the double-dot system. Fig. [5] shows that with increas- 
ing field, the double-dot conductance initially increases, 
then peaks at its unitary limit G = 2e 2 /h for a field 
value B** that for large A approaches 2Tk from above, 
and finally drops back toward zero for B 3> B** . The 
field B** is distinct from that characterizing the peak in 
7tA(0)Ai(0,0). In general B* < 2T K < B*\ but these 
three scales converge for A 3> A 2 . 



FIG. 6: (color online) Linear conductance G vs dot-2 level en- 
ergy £2 at zero temperature for Ui = — 2e\ = 0.5, A = 0.0627, 
and five different magnetic field values. The conductance is 
symmetric about the point £2 = of particle-hole symmetry. 
In nonzero fields, G peaks at some |ea| 7^ 0, apart from the 
special case B — B** ~ 0.045 ~ 2Tk, for which the conduc- 
tance is maximal at £2 = (as already seen in Fig. [SJ. Inset: 
Conductance vs magnetic field B for £2 = and 0.02. 



The initial rise in G with increasing field can be at- 
tributed to the progressive suppression of the Kondo ef- 
fect allowing dot 1 to become partially polarized and re- 
ducing the destructive interference between the Kondo 
resonance and the dot-2 resonant state. This change 
takes place — in agreement with the evolution seen in 
7rA(0) j4i(0, 0) and Mi — over a field scale that increases 
with A but is not just a constant multiple of Tk- By the 
point that the conductance reaches its peak at B = B** , 
the interchannel interference is clearly constructive since 
Eq. (|3"4"1) would predict a much lower conductance for dot 
2 alone. At still larger fields, the destruction of the Kondo 
resonance becomes complete and the dot-2 resonance is 
shifted far from the Fermi level, leading to a decrease of 
the conductance. 

Figure [5] illustrates aspects of the transport away from 
particle-hole symmetry. The main panel shows the vari- 
ation of the T = linear conductance at several differ- 
ent fixed magnetic fields as the value of £2 is swept by 
varying the voltage on a plunger gate near dot 2. For 
B = 0, the conductance increases from zero at £2 = 
and approaches 2e 2 /h for |£ 2 | 3> A 2 as the dot-2 res- 
onance is tuned away from the Fermi energy, thereby 
decreasing the destructive interference described in the 
previous paragraph and permitting perfect conduction 
through the Kondo many-body resonance^! For fixed 
B > 0, competition between Zeeman splitting of the 
dot-2 resonance and partial destruction of the Kondo ef- 
fect leads in most cases to an initial rise in G for small 
I £2 1 followed by a fall-off at larger \e%\. As the magnetic 
field increases from zero, the conductance peaks initially 
move to smaller |£ 2 |, then merge into a single peak at 
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G = 2e 2 /h for B = B** ~ 0.045 for the case A = 0.0627 
here, before separating and moving to larger \e 2 \ as B 
moves to still higher values. Thus B** can in principle 
be located as the only field at which G has a single peak 

VS £ 2 . 

The inset to Figure [6] compares the field variation of 
G(T = 0) for £ 2 = and for e 2 = 0.02. It is only in the 
former case (i.e., under conditions of strict particle- hole 
symmetry) , that the conductance has a single peak vs B 
and attains G = 2e 2 /h, whereas for e 2 ^ one finds a 
pair of peaks at G < 2e 2 /h. The presence of the unitary 
peak under field sweeps can therefore be used to identify 
the particle-hole-symmetric point in experiments. 

To better understand these results, we again turn to 
the noninteracting case U± — 0, where the linear con- 
ductance can be calculated by substituting the noninter- 
acting Green's function given by Eqs. ((29)) and ((301) for 
the full Green's function Q a in Eq. ([33]). At T = 0, this 
results in a conductance contribution 



G — *S G a 



p 2 



[l + e 2 a ][l + (e 2a -e lrT ) 2 } 



(35) 



where e\ a and e 2a are defined after Eq. (|32p . Equation 
([55)) correctly reduces to Eq. ([M)) in the limit \e\\ — > oo 
where dot 1 can play no role in the conductance. At 
particle-hole symmetry (ei = e 2 = 0), Eq. ([3"5)l gives 



G 



2e 2 



[S/2A(0)] 



h [1 + (B/2A 2 ) 2 ]{1 + [B/2A 2 



5/2A(0)] 2 } 

(36) 

which peaks at G = 2e 2 /h for B = B** = 2A, a char- 
acteristic field greater than the one B* — 2^/\ 2 — A 2 , 
at which 7rA(0) j4i(0, 0) reaches 1. Since we have seen 
above that G{T = 0) for the interacting case at particle- 
hole symmetry reaches 2e 2 /h for some B** > 2T K > B*, 
with B** — > 2T K for large A, the field dependence of the 
conductance reinforces the parallels between the large-A 
interacting problem and the noninteracting limit, with 
the many-body scale Tk playing the role of a renormal- 
ized A. 

An interesting feature of Eq. ([35)) is that it predicts 
conduction contributions Gf ^ G^ when particle-hole 
symmetry and time-reversal symmetry are both broken. 
In particular, for £\ > (or e\ < 0), the conductance 
polarization measured by 



G*t — Gi 
Gf + Gi 



(37) 



grows from ?/ = for B = to reach 77 = 1 (or 77 = —1) 
for B = 2|ei|, at which field e±i — (or £i-|- = 0), before 
decreasing toward zero for still larger fields. By contrast, 
keeping e\ = but allowing e 2 ^ results in variation 
of 77 with field, but does not allow one to achieve perfect 
polarization of the conductance. 

Spin-dependent conductance is also exhibited when dot 
1 has strong interactions. Figure [3a) shows the vari- 
ation of r\ with the dot 2 level energy e 2 in different 



T| o 




(b) 
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FIG. 7: (Color online) Conductance spin-polarization 77 (a) vs 
dot-2 level energy £2 at six fixed magnetic fields B, and (b) 
vs B for two values of £2. All data are for Ui = —2ei — 0.5, 
A = 0.0627, and zero temperature. In (a), r\ is odd about 
the point £2 = of particle-hole symmetry. Complete spin 
polarization of the conductance is achieved in the case B = 
0.01. Panel (b) shows a strong, nonuniversal variation of r\ 
with B for different values of £2. 



fields B 7^ for a symmetric dot 1 (XJ\ = ~2ei) and 
fixed A. The conductance spin-polarization is odd about 
the point e 2 = of particle-hole symmetry. For fields 
B < 2Tk — 0.042, 77 has the same sign as e 2l whereas for 
B > 2Tk, f] and e 2 have opposite signs. For each field 
value, |?7| peaks at a nonzero value of |£2|. One sees that 
a field B = 0.01 combines with a level energy \e 2 \ ~ 0.025 
to achieve complete destructive interference of the con- 
duction for one spin species, allowing passage only of a 
fully spin-polarized current through the device. The fact 
that reaching \r/\ = 1 in this manner — by varying e 2 while 
dot 1 is held at particle-hole symmetry (£1 = 0) — is im- 
possible in the noninteracting case U\ = indicates that 
the interference effects are more complex in the presence 
of strong interactions. 

Figure [7jb) shows the variation of r\ under field sweeps 
at two different values e 2 > 0. For each position of the 
dot-2 level, 77 changes sign at a nonzero B. For e 2 = 0.02, 
77 reaches +1 at a small field and then dips to nearly 
— 1 at a larger field before increasing back toward zero. 
For £2 = 0.1, by contrast, a small positive peak in 77 is 
followed at larger fields by a dip at (or very close to) — 1. 
This nonuniversal behavior reflects the subtlety of the 
interplay between the field and particle-hole asymmetry 
in controlling the constructive or destructive interference 
between transmission of electrons directly through dot 2 
and paths involving one or more detours to dot 1. 

Similar "spin-filtering" effects in a magnetic field have 
been investigated previously^ in the context of a single- 
mode wire, coupled near its midpoint via a tunnel junc- 
tion to a quantum dot (the "side dot"). In that work, 
conductance polarizations 77 = 1 and 77 = —1 (in the lan- 
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FIG. 8: (Color online) Linear conductance G vs dot-2 level 
position £2 for U\ = —2ei — 0.5 and A = 0.0627 at four 
temperatures T for (a) B = 0, and (b) B = 0.045 ~ 2T K0 . 
Temperatures are expressed as multiples of Tko = 0.021. 



FIG. 9: (Color online) (a) Conductance spin-polarization 77, 
and (b) conductance G vs dot-2 level energy £2 at different 
temperatures for Ui = -2ei = 0.5, A = 0.0627, and B = 0.01. 
Temperatures are expressed as multiples of Tko = 0.021. 



guage of the present paper) were found to occur at val- 
ues of the dot energy e d (r) = l) and £d(r/ = — 1) differing 
by a large scale exceeding the dot Coulomb interaction 
strength U. The combined effect of the applied magnetic 
field and the on-site interaction was a spin splitting so 
large that all traces of the Kondo effect were suppressed. 
These behaviors should be contrasted with those found 
here, where the £2 values that lead to 77 = ±1 differ only 
by an energy of order A2 (much smaller than U\). What 
is more, the complete spin filtering achieved in our setup 
depends crucially on the presence of Kondo many-body 
correlations. This point will become particularly clear in 
the next section, where we consider the effect of nonzero 
temperatures. 



B. Nonzero temperatures 

To this point, only zero-temperature results have been 
presented. This subsection addresses the effect of finite 
temperatures on the zero-bias conductance G and its spin 
polarization 77. Throughout the discussion, temperatures 
are expressed as multiples of a characteristic many-body 
scale Tko = 0.021, the system's Kondo temperature for 
£2 = 0, B = 0, and the representative value A = 0.0627 
that we have used in all our T > calculations. 

Figure M plots G vs £ 2 for A = 0.0627 in fields B = 
(panel a) and B ~ B** ~ 2T K0 (panel b). For B = 0, the 
effect of increasing temperature is a progressive suppres- 
sion of the Kondo effect and hence of the conductance 
channel involving the many-body Kondo resonance. As 
a result, G rises near £2=0 due to a lessening of the 
destructive interference between the Kondo channel and 
the single-particle resonance on dot 2 (discussed above 
in connection with Fig. [6]), but there is a decrease in the 
conductance at [£2! > A2, which is dominated by trans- 



mission through the Kondo channel. This trend results 
in a conductance peak at some |£2| 7^ for temperatures 
< T < Tko, which evolves into a peak centered at 
£2 = for T > Tko, m which regime transmission is 
dominated by the single-particle, Lorentzian-likc contri- 
bution from dot 2. 

Figure [8fb) reveals a very different behavior for B = 
B** ~ 2Tko- As described above, the T = conductance 
attains its unitary limit at £2 = due to constructive 
interference between the Kondo and single-particle con- 
ductance channels, and G decreases monotonically with 
increasing |e 2 |. Raising the temperature over the range 
T < Tko leads to suppression of the Kondo conductance 
channel but has little effect on the single-particle channel, 
leading to a decrease in G that is strongest for £2 = 0. 
Once the temperature passes Tko, the variation of G with 
£2 increasingly reflects the field splitting of the dot-2 en- 
ergy level, with peaks centered at £2 ~ ±^B. 

The influence of temperature on the spin polarization 
of the conductance is shown in Fig. [Hta), which focuses 
on the case B = 0.01 that we know from Fig. [7] yields 
full spin polarization (7/ = ±1) at zero temperature for 
£2 — ±0.025. As T increases from zero, the peak spin po- 
larization is lowered, presumably due to a combination 
of two effects: (i) a reduction in the destructive interfer- 
ence between the Kondo and single-particle conduction 
channels for one spin species a leading to an increase 
in — Im7^-(w = 0,T) entering Eq. ([3]); and (ii) thermal 
broadening of —df/doj in Eq. (|5|) leading to sampling of 
w values having nonzero -Im7^j(w,T = 0). At higher 
temperatures, T ~ Tko, the suppression of the Kondo 
conductance channel unmasks oscillations in 77 vs £2 that 
result from shifts in the spin-resolved energy levels in 
dot 2. These oscillations are much less pronounced than 
the polarization variations at lower temperatures and the 
maximum values of 77 1 are about an order of magnitude 
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FIG. 10: (Color online) Effect of a nonzero dot-2 interaction 
(U2 > 0) on (a) the conductance spin-polarization 77, and (b) 
the conductance G, both plotted vs dot-2 level energy 62 for 
the same parameters as in Fig. [9] Open symbols correspond 
to U2 — and filled symbols to U2 = 0.01. Temperatures are 
expressed as multiples of Tko = 0.021. 



smaller than those obtained in the Kondo regime. The 
above trends with increasing temperature serve to un- 
derline the central role of many-body effects in achieving 
complete or near-complete spin filtering, in stark contrast 
to the situation considered in Ref. l32l . 

Figure EJb) plots the total conductance G vs £2 corre- 
sponding to each of the r\ vs £2 traces in Fig. EJa). There 
is a close correlation (although not a perfect match) be- 
tween the £2 values of the peaks in G and of those in 
77 1 . This suggests that measurements of the total con- 
ductance can provide a useful starting point for experi- 
ments seeking to optimize the system's spin-filtering per- 
formance. 

In this paper, we have focused on the special case 
U2 = that facilitates analysis through mapping of the 
double-dot problem to a single-impurity Anderson model 
with an energy-dependent hybridization function. How- 
ever, the conductance features described above by no 
means depend on the condition U2 = 0. In fact, qual- 
itatively similar results are obtained for an interacting 
dot 2 provided that U2 is small compared to the level 
broadening A2. This is illustrated in Fig. [TUl which 
compares the £2 dependence of the conductance and of 
its spin-polarization for U2 = (data from Fig. [9} and 
U 2 = 5A2 = 0.01, both for the lowest (T = 0) and high- 
est (T = QA60Tko) temperatures shown in Fig. [9] Apart 
from a small shift in the point of particle-hole symme- 



try, which moves from £2 = to £2 = —\U2^ the other 
essential features (such as the complete spin polarization 
at zero temperature) are unaffected by the presence of 
Coulomb repulsion within dot 2. 

VI. CONCLUSIONS 

In this work, we have investigated the effect of an ap- 
plied magnetic field on a strongly interacting "hanging 
dot" side-coupled to external leads via a weakly interact- 
ing "bridge dot." Our numerical renormalization-group 
results show that the interplay of electronic interference, 
the Kondo effect, and Zeeman splitting brings about 
qualitative changes in the spectral and transport proper- 
ties of this system compared to those in the well-studied 
case of a single quantum dot. 

We have found, for instance, that the value of the in- 
teracting dot's zero-temperature spectral function at the 
Fermi energy does not decay monotonically with increas- 
ing field. Instead, the presence of an extra energy scale 
related to the interdot coupling introduces nonuniversal 
behavior, and in some cases leads to the appearance of 
one or two maxima in the Fermi-energy spectral func- 
tion at nonzero values of B. These features can be ex- 
plained by the presence of a parameter-dependent phase 
in the Fermi-liquid relation expressed by the Friedel sum 
rule, generalized to account for an energy-dependent hy- 
bridization function. 

One of the signatures of the interplay of site and 
spin degrees of freedom in this double-dot device is the 
appearance of spin-polarized currents between the two 
leads. We have shown that the degree of spin polariza- 
tion can be tuned up to 100% by changing gate voltages 
and/or small magnetic fields in the system. These results 
underscore the flexibility of quantum-dot systems for ex- 
ploration of novel effects in correlated electron physics. 
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